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ABSTRACT 

We describe the results of a search for microlensing events affecting stars 
in the outer bulge and inner disk of M31, due both to masses in M31 and the 
Galaxy. These observations, from 1994 and 1995 on the Vatican Advanced 
Technology Telescope and KPNO 4m, are sufficient to rule out masses in the 
range of ~ 0.003 Mq to 0.08 Mq as the primary consistuents of the mass of 
M31 towards this field. Furthermore we find six candidate events consistent 
with microlensing due to masses of about 1 Mq, but we suspect that some of 
these may be cases where long-period red supergiant variables may be mistaken 
for microlensing events. Coverage from anticipated data should be helpful in 
determining if these sources maintain a constant baseline, and therefore are best 
described by microlensing events. 

Subject headings: galaxies:individual(M31) - dark matter - gravitational lensing 
- star:variables:other 
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1. INTRODUCTION 

One of the most significant and stubborn mysteries in astrophysics today concerns 
the nature of the dark matter in spiral galaxies (e.g. Rubin et al. 1978). The least radical 
candidate for such dark matter is baryonic objects which are too large to be detected as 
dust or gas. Indeed, if the Hubble constant is not too large, a significant fraction of the 
baryons, as implied by Big Bang nucleosynthesis, must be hidden as dark matter (Walker 
et al. 1991). Arguments have been made for why such dark baryonic objects are unlikely 
on individual mass scales of atoms to brown dwarfs (Hills 1986, Hegyi & Olive 1986). Still, 
objects of primordial composition and more massive than about lO^^M© might be expected 
to resist evaporation until the present day (de Riijula et al. 1992), while masses smaller 
than about 0.077 Mq would fail to ignite as stars (Burrows et al. 1993). 

Astrophysicists' frustration explaining the dark matter with any directly detectable 
objects has led to the suggestion that gravitational microlensing might be used to at least 
betray the presence of individual objects via their effects on background stars as sources 
(Paczynski 1986), and thereby give some indication as to their mass. Such searches have 
recently taken place towards the Large Magellanic Cloud (LMC) (Alcock ct al. 1996, 
Aubourg et al. 1995) and Bulge (Alcock et al. 1995, Paczynski et al. 1994), with searches 
towards the LMC ruling out most of the dark matter being composed of substellar-mass 
objects (Aubourg et al. 1995, Alcock et al. 1996) heavier than about 10~^Mq, while 
suggesting that a large fraction might have the same component mass as low-mass stars 
(Alcock et al. 1996). Given the uncertainty of the Galactic halo's distribution of MAssive, 
Compact, Halo Objects (MACHOs) and therefore the lensing geometry leading to events, 
the relationship between mass and observed microlensing lightcurve timescale is still 
unclear. 

In part because of its unique geometry with respect to Earth and partially due to high 
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predicted optical depths (r) due to lensing, M31 is a uniquely powerful venue for studying 
microlensing. Early we realized that an M31 microlensing survey would show particular 
advantages if the practical aspects of studying such a distant, crowded field of stars could 
be overcome. 

We found such an approach, briefly outlined by Crotts (1992) with a complete 
description of the realistic technique and preliminary results found in Tomaney and Crotts 
(1996, hereafter TC). By subtracting images in a time sequence, then performing "difference 
image photometry" (DIP, also know as "pixel lensing"), we can study the residual point 
sources due to variables, while the signals from the many crowded, non-varying stars 
subtract away. With a practical method of observation and analysis, we can exploit the 
advantages inherent in studying M31: 1) very small component mass limits, due to the 
small angle subtended by the photosphere of M31 stars compared to the Einstein radius of 
objects of solar mass (c.f. TC for low- mass results), 2) the abihty to study different parts 
of M31, thereby studying the spatial distribution of microlensing objects, 3) the abihty to 
study many stars at once in fields of high r, thereby detecting events in short periods of 
observation, and 4) the constrained microlensing geometry, due to the fact that lensing 
mass is concentrated over the center of the galaxy, thereby allowing a better determination 
of the MACHO mass given microlensing event timescale. 

It is our hope that by studying M31 in this way, both its halo and bulge, one can 
more readily understand both these results and those obtained in the Galaxy. This paper 
presents our results from our first season of observation toward this goal. 
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2. Observations 

Observations were made primarily at the Vatican Advanced Technology Telescope on 
31 December 1994, 31 January, 17 October through 4 November and 20 November through 
December 3, 1995, using an imager especially made for accommodating DIP. Additional 
observations were made with the prime focus KPNO 4-m Prime Focus Camera on 24-27 
September 1994 and 28 August 1995. 

The VATT data were taken using field centers of 00:43:16.5 +41:11:33 and 00:42:13.4 
+41:20:44 (J2000), and rotation of the CCD so that its sides were roughly parallel to M31's 
principle axes. The first field looks past the bulge, intercepting the far side of the disk along 
the minor axis. Most data were obtained in the first field, such that light curves in the 
second are too incomplete to treat here. Observations on the KPNO 4m were centered close 
to the first VATT location, such that the 4m field (16.4 arcmin on a side) encompasses all 
but 0.1% of the VATT field (11.3 arcmin on a side). The two telescopes' fields are rotated 
38° with respect to each other. We do not consider here 4m data faUing outside the VATT 
field. Assuming a distance to M31 of 770 kpc (4.46 arcmin/kpc), the VATT field covers a 
range of 0.4 to 2.9 kpc along the minor axis, which projects to 1.8 to 13.1 kpc along the 
disk, assuming an inclination i — 77°. We use the filter bands described in TC, essentially 
broad R and I bands to match the fact that most of our target stars in M31 are red giants 
or super giants and therefore brighter in redder bands. 

3. Analysis 

The process of difference image photometry consists of careful flat-fielding, coordinate 
registration and photometric scaling of the data, followed by point spread function (PSF) 
matching between frames (detailed in TC). Using a convolution kernel approximating the 



-6- 



quotient between PSFs in the image from a particular epoch and that from a high-5'/A^, 
good-seeing stack of many images, the PSFs are matched either by degrading the image 
stack or single epoch's image to match the other. In the case of the VATT data, the entire 
frame can be corrected in this way due to the optical design of our imager, which insures 
that a single convolution kernel is need over the entire (or nearly all) of the image (see TC 
for more details). This is especially easy because in these data we deal only with nightly 
sums, composed of ~ 5 — 20 similar exposures, so that individual irregularities of single 
exposures, such as guiding errors, average out. After PSF matching is accomphshed, the 
stack (taken from many night's data) and the individual night's image are subtracted, 
leaving a field of noise at nearly the photon shot noise level, as well as isolated positive 
or negative point source residuals due to variable stars. An example of this is shown in 
Figure 1, for the fourth candidate microlensing event detailed in Figure 2. DIP is completed 
by performing aperture photometry on these isolated residual sources, which can then be 
incorporated into light curves. 

The error bars presented with the lightcurves in Figures 2 and 3 show the fluctuations 
in the difference image on the scale of the PSF in regions adjacent to each residual source. 
Sources were catalogued by requiring at least a Aa detection in at least two nightly sums 
(or 6(7 in the 24% of the image containing the bulge and closest to the minor axis), then 
tracing the lightcurve by aperture photometry in other epochs at the same location. In 
future papers we will track additional sources by 1) sampling variations on sub-night scales, 
and 2) summing difference images so that weaker residuals can be tracked over longer 
timescales. Even without these refinements, however, we locate over 2000 sources within 
the VATT field. 
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4. Results 

The results of this construction of hghtcurves from nightly sums is that no source is 
only on two consecutive nights, and that none of the sources seen, with one exception, is 
consistent with microlensing events on any but nearly the longest timescales sampled by 
our survey. In the latter cases, we portray the lightcurves of the six candidate events in 
Figures 2, and other information in Table 1, including their positions (J2000) and distance 
along M31's minor axis (d). Assuming that they arc microlensing events, several other 
parameters can also be extracted: the duration (Einstein radius crossing time te), lensing 
impact parameter (normalized to the Einstein radius: Uo — u/Re), and source baseline 
magnitude (R). Not given are the two other fit parameters, time of peak amplification and 
flux zero-point offset due to image subtraction. Additionally we give the goodness of the 
best lensing fit (for point sources and masses), and the most probable mass of the lens. It 
appears that fit residuals are slightly larger than expected from photometric measurement 
error alone, seen particularly as a surplus in the number of 3o" or greater residuals, which 
are inconsistent with neighboring points. One possibility for this noise is underlying RR 
Lyrae variables, which should be evident at the 1-2(7 level, either coincident with the source 
or in its photometric background annulus. We will investigate this problem further in 
Tomaney et al. (1996). We stress that we do not claim that these are microlensing events at 
least until their lightcurves are observed to fall and remain at the pre-event baseline during 
the 1996 observing season or thereafter. 

One reason for our caution is the lightcurve shown in Figure 3 for the variable star 
found at 0M3"*37.^9 +41°14'57" (J2000), 1.92 kpc from the major axis. It is fit well 
by a microlensing lightcurve during its rise and fall, but does not maintain a consistent 
baseline before and after the event. Upon inspection of lightcurves of Mira-type variables 
(Wesselink 1987), we find a small fraction whose lightcurves around maximum light mimic 
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the behaviour of microlensing lightcurves. A sparse sampling of points beyond maximum 
light, if chosen unfortuitously, might fail to distinguish such a variable from a microlensing 
event. We find several such variables in our VATT field. Further reasons for suspicion is 
the similarity in timescales to those of miras (except the first event), and similar shapes, 
indicated by Ug values which cluster around 0.6 (except for the second event). Additionally, 
it is strange that all sources have R ^ 21, close to the magnitude that would correspond to 
a mira pulsation (given the inferred Uo), but brighter than what we might expect for lensed 
sources given the luminosity function of stars in the field. We suspect that several of these 
events are not due to microlensing at all, but might be associated with bright variable stars. 
Another season of observation, which we plan, will determine if these sources maintain a 
constant baseline and are therefore hkely to be lensed. 

The reality of these events can be tested in terms of the distribution of Ug values via a 
one-sample Kolmogorov-Smirnov (K-S) test. The theoretical distribution is derived using a 
luminosity function in i? oc 10"^, where a — 0.59 best describes the behavior of star counts 
and surface brightness fluctuations in the field, and agrees with other works (Tomaney & 
Crotts 1996, and references therein). Very low Ug values (mo<0.02) are not realistic since 
high amplification events have timescales too short for us to detect; high Ug events {ug^3.5) 
have amplifications too subtle for us to detect, as well. We take as the Ug upper bound the 
lesser of the above upper bound and the maximum value providing sufficient amplification 
to reach our flux threshold for a given magnitude. For values of Ug between these limits, we 
assume a uniform distribution of events in Ug, at a, given magnitude. The largest value of 
the K-S distance D occurs at the smallest observed Ug = 0.369, due to the lack of small Ug 
events, and has a value D ^ 0.7. Assuming that all six candidates are true microlensing 
events, the null hypothesis (consistency with microlensing) is rejected at the 99.5% level. If 
half the candidates are microlensing events (and the minimum still Ug = 0.369), the null 
hypothesis is rejected at approximately the 90% level. It is unlikely that all of the events 
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are due to microlensing, but this test cannot rule out that a large fraction may be. 

These caveats aside, if these events are microlensing events, then we can say several 
things about them. The first and second events land in the bulge-dominated region, and 
hence likely involve bulge sources. The third and fourth might be due to disk sources (but 
have a high probability of belonging to the bulge) and also rest in the region where bulge 
lenses may dominate over halo lenses. The fifth and sixth events, if genuine, might easily 
be halo lenses acting on disk sources. In the case of third through sixth events, the most 
probable source-lens distance is d/cos i, allowing us to compute a most likely mass, given 
te- (We assume a disk rotation speed of 260 km s~^ [Braun 1991] and a halo/bulge velocity 
dispersion of 160 km s"^ [Kent 1989], of which ^2/3 is in the transverse direction. Earth's 
transverse motion is negligible.) We assume no rotation of the bulge; it could be as large 
as ~100 km s~^ in our field (Kent 1989), meaning that inferred masses might tend to split 
into a bimodal distribution of under- and overestimated values, with peaks differing by as 
much as a factor of two in timescale, or four in mass. 



5. Discussion 

Several approaches have been taken to estimating the predicted r in M31 due to its 
own mass distribution. Initially Crotts (1992) just approximated the entire mass of M31 
as an r"^ density distribution, which produces an optical depth for far-side disk stars of 
r ^ 10"^. The presence of a core saturation radius will reduce r in the center of M31 
while maintaining this high plateau value at larger radii. Jetzer (1994) considers the 
effects of only the dark matter halo, with a large core radius of 5 kpc, and finds a value of 
T — Ix 10~^ in the center of M31, rising to 3 x 10~^ at the outer edge of our field. Han and 
Gould (1996) treat both the halo and bulge of M31 and find r = 7 x 10~^ in the center, 
dominated by the bulge, falling to 3 x 10~^ at the outside edge of our field, where the halo 
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dominates. There is a factor of about 1.3 disagreement between the lower values of Jetzer 
versus Crotts and Han & Gould due to different assumed values of M31's rotation velocity. 
Lensing of disk stars by other disk stars produces a r component of 4 x 10~^ (Gould 1994), 
while a standard Galactic halo model adds r f« 1 x 10"^ (Paczynski 1986). Together, these 
components sum to at least r pa 5 x 10~^ throughout the field, which is the value that we 
will adopt for the sake of discussion. Note that this is about an order of magnitude greater 
than that suggested by Galactic survey results towards the LMC (to6s = 2.9lo!9 x 10~^, 
Alcock et al. 1996), or for predicted Galactic/LMC halo values {Tmodei — 4.7 x 10~^, Alcock 
et al. 1996), but only slightly larger than Galactic Bulge results [Tobs = (3.3 ± 1.2) x 10~®, 
Paczynski et al. 1994; Tobs = (3.9 ± 1.8) x 10"^ Alcock et al. 1995). 

Prom our previous constraints on the luminosity function of stars in our field (TC), we 
have estimate that we are sensitive to detectable microlensing of any of 6.9 x 10^ stars in our 
field. These data are primarily sensitive to timescales ranging from 2^ to 10*^, corresponding 
to 0.003 Mq to 0.08 Mq. We have 13 and 2 sample times respectively corresponding to 
9.0 X 10^ and 1.4 x 10® star-epochs. The predicted number of events for this mass range 
given a Toai+Msi of 5 x 10~® is 45 to 7 events. Except for one possible detection at the upper 
end of this range, we find no events on these timescales, thereby eliminating this mass range 
as a 100% contribution to the mass of M31 at considerably better than 95% confidence. On 
the other hand, we expect to detect approximately 2 events (given 100% efficiency) if the 
mass of M31 is made entirely of 1 Mq objects, while we see six candidates, half of which are 
at this scale or larger. This argues that some of these may not be caused by microlensing. 

Our primary result is (1) the lack of any detection corresponding to masses up to 
0.08 Mq (with perhaps one exception), and (2) the possible detection of events on the scale 
of about 1 Mq. The number of such events on this larger timescale, however, is significantly 
greater than would be predicted given models of the lensing optical depth, so might indicate 
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contamination by variable stars. Both of these results are consistent with microlensing 
searches in both the Bulge and halo of our Galaxy (e.g. Paczynski et al. 1994, Alcock 
et al. 1995, 1996) in which few, if any, substellar masses are detected. Likewise, slightly 
sub-solar masses are indicated as the primary cause of microlensing events both towards the 
LMC and Bulge. A further season of data will determine whether our six candidates are 
simply variable stars, or exhibit constant baseline, implicating them as microlensing events. 

We appreciate the assistance of Robert Uglesich in the preparation of this paper. We 
also thank the Vatican Observatory Research Group, in particular Richard Boyle and Chris 
Corbally, in assisting our efforts on the VATT, and the Steward Observatory engineering 
group, in particular Richard Cromwell and Steve West. Ata Sarajedini aided our program 
at the KPNO 4m. This research was funded by the David and Lucile Packard Foundation. 
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Figure 1: The left panel shows reduced but unsubtracted subimage (about 450 arcsec^) prior 
to DIP. (Actually, the average of our 24 VATT epochs from autumn 1995 is shown.) Some 
of the sources shown here are variable, but at a level far below the average surface 
brightness fluctuation in the image and below a level which can be distinguished by eye. 
The right panel shows the difference image of the same region, after a high-5'/A^ sum of 
images has been scaled, registered, PSF adjusted and subtracted using DIP, for UT 17 Oct 
1995. Several variable sources are revealed by the image subtraction, some fainter than 
the average (black), some brighter (white). The circled source corresponds to the fourth 
microlensing candidate shown in Figure 2. 

Figure 2: the hghtcurves of the six candidate microlensing events described in the text and 
Table 1. One count is equivalent to i? = 31.09. 

Figure 3: the example described in the text of one of several hghtcurves well-flt in the peak 
by a microlensing model, but which does not maintain a consistent baseline outside of the 
peak. 
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Table 1. Candidate M31 Microlensing Events 



RA 


Dec 


Minor Axis 


Duration, 


Probable 


Impact 


Baseline 




(J2000) 


(J2000) 


Distance, 


te (days) 


Mass, 


Parameter 


Magn., 








d (kpc) 




m (Mq) 


vs. Re, Uo 


R 






+41°14'27" 


0.59 


7.8 


(0.09)" 


0.648 


21.08 


0.75 


0M2"*42^3 


+41°11' 2" 


0.62 


55.7 


(4.3)« 


0.369 


21.14 


1.67 


0M2"^54M 


+41°10'55" 


1.02 


32.6 


0.90'' 


0.680 


20.22 


1.35 


0M3™14^8 


+41°12'32" 


1.49 


39.3 


0.90^ 


0.501 


20.77 


1.62 




+41° 5'52" 


2.67 


31.7 


0.32 


0.590 


20.93 


2.23 


0M3"'49^0 


+41°11'28" 


2.77 


39.9 


0.49 


0.690 


20.61 


2.01 



^Source probably in bulge, so mass estimate is unreliable. 
"^Source in bulge or disk; mass estimate assumes disk. 



